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Narrative  Description  of  the  jarogram: 


The  goal  of  this  program  is  to  develop  a  basis  for  developing  advanced  crack  growth  and  service  life 
prediction  technologies  for  predicting  the  service  life  of  solid  rocket  motors.  The  objectives  of  this 
program  are  to  (1)  gain  a  fundamental  understanding  of  fracture  and  crack  growth  behavior  in  solid 
rocket  motors;  (2)  investigate  the  effects  of  damage,  material  nonlinearity,  pressure,  and  loading  rate 
on  crack  growth  behavior  in  a  solid  propellant;  (3)  simulate  crack  growth  behavior  and  gain  insight 
for  improving  crack  growth  resistance  in  solid  propellant;  and  (4)  determine  the  crack  growth 
behavior  of  an  interfacial  crack  in  bimaterial  bond  systems.  The  main  issues  in  service  life  prediction 
of  solid  rocket  motors  are  the  lack  of  a  fundamental  understanding  of  crack  growth  behavior  under 
service  loading  conditions  and  a  reliable  methodology  to  predict  crack  growth.  The  main  technical 
challenges  are  microstructure  effects  on  damage  initiation  and  evolution,  large  and  tim^dependent  -f 
deformation,  short  crack  and  stress  raiser  interaction,  and  multi-layer  structures  with  time-dependent 
material  properties  and  property  gradients.  The  program’s  basic  approach  involves  a  blend  of 
analytical  and  experimental  studies.  In  general,  mechanisms  and  mechanics  involved  in  cohesive 
fracture  in  a  solid  propellant  and  adhesive  fracture  in  bond  systems  are  emphasized.  In  this  program, 
nonlinear  viscoelasticity,  fracture  mechanics,  experimental  mechanics,  damage  mechanics, 
nondestructive  testing  and  evaluation,  and  numerical  modeling  techniques  will  be  used. 

These  research  studies  address  a  number  of  important  subjects  such  as  cumulative  damage  and  crack 
growth  behavior  in  solid  propellants,  statistical  nature  of  crack  growth,  and  bonded  interface  failure. 

The  results  of  these  studies  have  the  potential  of  becoming  some  of  the  most  significant  contributions 
to  the  rocket  industry  and  research  community. 


Detailed  Technical  Approach  for  Next  Fiscal  Year: 

In  FY  01  there  are  four  major  tasks:  Task  1-  investigating  the  effects  of  pressure,  strain  rate,  and 
damage  on  short  crack  growth  behavior  in  a  solid  propellant,  Task  2  -  determining  the  validity  of 
using  linear  fracture  mechanics  to  characterize  the  crack  growth  behavior  in  a  solid  propellant,  Task  3 
-  including  pressure  effect  in  numerical  modeling  of  crack  growth,  and  Task  4  -  determining  the  strain 
rate  effects  on  the  constitutive  behavior  as  well  as  on  the  local  damage  and  deformation  near  the  tip  of 
interfacial  cracks  of  a  propellant/liner/propellant  bond  system. 

\  tH  / Task  1 :  Investigating  the  Effects  of  Pressure,  Strain  Rate,  and  Damage  on  Crack  jrowth  Behavior  in  a 
*  Solid  Propellant  ~ 


<V" 

In  Task  1,  experiments  were  conducted  on  uniaxial  specimen  jwith  two  different  initial  crack  lengths 
(0.1  in.  and  0.3^in.)  at  four  different  displacement  rates  (0.2rni/min?  2  in/min,  50  in/min,  and  200 
pis  in/min)  under  three  confined  pressures  (ambient,  500  psi,  and  1000  psi).  Preliminary  results  of  data 

+o  ^  —  analysis  revealed  that  cracks  grew  under  confined  pressure.  The  decrease  in  crack  growth  rate 
under  confined  pressure  is  due  to  the  suppression  of  the  development  of  damage  in  the  material  and 
the  decrease  in  strain  gradient  near  the  crack  tip  (Fig^^^. the  time  of  this  writing,  the  basic  damage 
vl  mechanisms,  including  damage^ initiation  and  evolution,  uncllr  multiaxial  (constant  strain  rate  and 

"  confined  pressure  conditional^ not  clear.  In  order  to  provide  a  fundamental  understanding  of  the 
juA X  fis  —  effects  of  multiaxial  loadings  on  crack  growth  behavioj^dditional  studie^including  experiments  and 
ft^^v^unierical  modeling^will  be  conducted  in  FY  02.  ^ 

^  During  FY  01,  in  addition  to  experimental  study  on  the  effect  of  confined  pressure  on  crack  growth 

y*  behavior,  research  was  centered  on  the  development  of  a  rate-independent  constitutive  model, 

including  the  pressure  effect.  Particularly,  the  study  was  conducted  on  modeling  the  effect  of 
hydrostatic  pressure  on  the  progressive  damage.  In  order  to  model  the  effect  of  hydrostatic  pressure 
on  the  damage,  the  damage  theory  previously  developed  was  modified,  and  it  was  included  in  the 
constitutive  model.  Also,  the  technique  developed  previously,  called  the  multi-scale  approach,  was 
used  along  with  the  damage  theory  at  the  matrix  material  level.  The  multi-scale  approach  is  based  on 
the  interconnection  between  the  micro  level  (i.e^particle  and  matrix  material  level)  and  the  macro 
level  (i.eyythe  particulate  composite  level),  and  it  was  used  to  predict  the  stress-strain  behavior  of  the 
material  under  ambient  pressure  and  1000  psi  confined  pressure  conditions.  An  excellent  agreement 
exists  between  the  predicted  and  the  experimental  stress-strain  curves  (Fig^2). 

Task  2:  Crack  Instability  and  Growth  Models 

Sub-Task  1:  Instability  Criteria  for  Short  Crack  Growth 

In  this  task,  the  instability  criteria  for  the  onset  of  growth  of  short  cracks  willbe  developed.  Based  on 
last  year’s  study,  we  found  that  for  short  cracks,  defined  as  the  crack  length j&^qual^or  less  than  0.1 
m.,  classic  fracture  mechanics  cannot  be  used  to  determine  the  critical  condition  for  the  onset  of  crack 
growth.  In  FY  01,  a  linear  fracture  mechanics  solution  was  modified  to  predict  the  behavior  of  short 
cracks.  An  effective  crack  length  was  introduced  into  the  Mode  I  stress  intensity  factor.  By  using  the 
effective  crack  length,  a  reasonably  good  agreement  exists  between  the  fracture  tougluiessfiSTor  the 
onset  of  growth  of  short  and  long  cracks.  Additional  work  will  be  conducted  in  F^Q2  to  develop  the ' 
instability  criterion  for  the  onset  of  growth  of  short  cracks,  based  *-n  continujum  and  fracture 
mechanics  theories. 
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Sub-Task  2:  Nonlinear  Crack  Growth  Model 


In  this  task,  nonlinear  and  linear  fracture  mechanics  together  with  numerical  modeling  techniques  will 
be  used  to  determine  the  capability  of  using  linear  fracture  mechanics  to  characterize  the  crack  growth 
behavior  in  a  solid  propellant  under  confined  pressure.  According  to  the  previous  studies,  under 
ambient  pressure,  the  crack  grew  rapidly  when  the  specimen  was  deformed  beyond  the  point  of 
maximum  tensile  stress.  However,  based  on  recent  experimental  observation,  it  was  noted  that,  under 
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onfined  pressure,  a  considerable  amount  of  stable  crack  growth  took  place  when  the  specimen  was 
deformed  beyond  jthe  maximum  applied  stress.  This  type  of  crack  growth  behavior  is  significantly 
different  from  thaMinder  ambient  pressure.  Therefore,  in  order  to  predict  the  criticality  of  a  crack  in  a 
solid  rocket  motbr  under  a  pressurization  condition,  a  nonlinear  fracture  mechanics  theory  is  needed 
to  characterize  the  crack  growth  behavior  and  to  develop  a  crack  growth  model.  In  FY  02,  efforts  will 
be  centered  on  nonlinear  analysis  of  crack  growth  behavior  and  determining  the  applicability  of  using 
linear  fracture  mechanics  to  characterize  the  crack  growth  behavior  under  confined  pressure. 

fk  ^  kj/is 

Sub-Taslc3:  Determining  the  Representative  Volume  for  the  Validity  of  Homogeneous  Continuum 
Assumption  of  the  Materials  -  Mtk  m 


In  FY  00,  real-time  x-ray  techniques  were  used  to  investigate  the  validity  of  the  homogeneous 
continuum  assumption  of  the  solid  propellant.  The  real-time  x-ray  test  data  indicated  that  the  particles 
were  randomly  distributed  in  the  particulate  composite  material.  When  the  size  of  an  area  becomes 
larger  than  2  mm  x  2  mm,  the  microstructure  becomes  statistically  uniform.  In  other  words,  the  area 
size  of  2  mm  x  2  mm  is  about  the  smallest  to  represent  the  average  smeared  behavior  of  the 
.  ^^composite.  This  2  mm  x  2  mm  area  is  defined  as  a  representative  area.  In  FY  01,  Lockheed-Martin 
~  Research  Laboratory’s  Computed  Tomography  system  was  used  to  determine  the  representative 
volume  of  the  material.  Preliminary  data  showed  that  the  representative  volume  of  the  material  is  2 
JUltks*  nm  x  ^  rnm  x  2  mm/" (Fig.  3).  This  information  can  be  used  to  correlate  the  material  property  to 
r*  f^noo(  nondestructive  testing  parameters,  providing  a  basis  for  conducting  stress  analysis  of  structures  from 
nondestructive  testing  data.  In  FY  02,  the  Computed  Tomography  System  will  be  used  to  monitor 
damage  initiation  and  evolution  processes  inside  the  specimen  under  incremental  strain  conditions 
and  to  correlate  microstructure  to  surface  strain  measured  on  the  specimen  surface. 


Sub-Task  4:  Investigating  the  Effect  of  Strain  Rate  on  the  Critical  Initial  C  rack  Length  in  the  Material 
West* 

In  FT^IO,  a  technique  was  developed,  based  on  fracture  mechanics  and  probabilistic  mechanics,  to 

predict  the  critical  initial  crack  size,  which  is  responsible  for  the  fracture  of  the  specimen.  In  FY  01, 

the  developed  technique  was  modified  to  include  rate  effect.  The  modified  technique  can  be  used  to 

predict  the  initial  critical  crack  size  at  different  strain  rates  with  good  accuracy.  For  example,  a^fk  v 

18.182  in/in/min  strain  rate,  the  predicted  and  the  measured  initial  critical  crack  siz^are  0.146  irg^nd  s 

,l>  j  —  0.12  ii^  respectively.  The  results  also  show  that  the  critical  initial  crack  size  is  insensitive  to  strain 

,  L_  rate.  For  example,  at  0.067  in/in/min  and  0.727  in/in/min  strain  rates  the  predicteddnitial  critical  crack 

10  L  _ ,  —  siz^are  0.13  i%and  0.1 19  iij,  respectively.  The  determination  of  the  initial  critical  crack  size  and  its 

statistical  distribution  function  will  make  statistical  and  reliability  analyses  of  crack  growth  feasible.  er  wj- 
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Task ! 


Numerical  Modeling  of  Time-Dependent  Crack  Growth  Behavior/-"  Jddc  "I*-  ja*-^ 


~i  this  task,  a  nonlinear  viscoelastic  constitutive  model  was  incorporated  into  a  finite  element  code 
(FEAP),  which  was  used  to  model  the  crack  growth  behavior  in  a  solid  propellant.  The  response  of 
the  undamaged  material  is  assumed  to  be  governed  by  a  homogenized  linear  viscoelastic  constitutive 
relation.  The  bulk  and  the  shear  moduli  of  the  material  are  assumed  to  undergo  degradation  with 
accumulation  of  damage.  The  critical  damage  parameter  is  taken  to  be  the  maximum  volume 
dilatation  accumulated  in  the  material.  In  the  finite  element  analysis,  the  crack  propagation  is 
simulated  by  using  a  node  release  technique  after  having  attained  a  critical  level  of  dilatation  in  the 


elements.  In  F^f)l,  the  effect  of  pressure  on  the  constitutive  behavior  was  included  in  the  constitutive 
model,  which  was  incorporated  in  the  FEAP  computer  program.  In  FY  02,  the  accuracy  of  the 
computer  program  will  be  determined  by  comparing  the  numerically  predicted  and  experimentally 
measured  crack  §rowth  behavior. 
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Task  4:  Interfacial  Fracture  of  Bimaterial  Bond  Systems 

"1W5  «  uLa  Ca.P  A  pwifclj 

.Sub-Task*  l^hotoelasticAnalysis  of  Three-Dimensional  Effects  of  Cracking  of  Motor  Grain 
Geometries  uimer^fniernalrressure  Jpads/--  ?U  aJL" 

'JcLJn-  fcr/hJt 

In  FY  01,  a  series  of  experiments  on  photoelastic  scale  models  of  motor  grain  were  conducted  using 
the  frozen  stress  method.  The  results  of  preliminary  data  analysis  revealed  that  symmetric  cracks 
penetrated  to  the  outer  surface  of  the  model  before  the  off-axis  cracks  had  grown  significantly.  Data 
also  revealed  a  shear  mode  along  with  Mode  I  in  the  off-axis  cracks  before  they  were  completely 
turned  (Fig.^T  In  F^2,  a  detailed  analysis  will  be  conducted  to  determine  the  three-dimensional 
effects  on  the  stress  distribution  at  various  locations  along  the  crack  front. 
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Sub-Task  2:  Deformation  and  Failure  Mechanism  of  Propellant/Liner/Propellant  Bonded  Specimens 


^  ^  -  In  FVt)l?X series  pf  experiments  on  propellant/liner/propellant  bonded  specimens  were  conducted  at 
0.01  in/min  displacement  rate.  A  computer  aided  speckle  interferometry  technique  was  used  to 
4+  determine  the  displacement  and  strain  distributions  in  the  specimen.  Two  interface  debonding  modes, 

debonding  from  the  center  and  the  comer  of  the  interface  of  the  specimen,  were  observed.  These 
debonding  modes  appeared  to  be  related  to  the  specimen  geometry.  In  addition,  the  strain  rate  in  the 
interphase  and  liner  layers  increase  with  increasing  time  and  are  significantly  different  from  the 
constant  applied  strain  rate  (Fig.5).  In  FY  02,  the  effect  of  the  applied  strain  rate  and  interfacial 
strength  on  the  deformation  and  strain  distributions  in  the  bonded  specimen  will  be  investigated. 

sub-Task  3:  Numerical  modeling  of  the  Propellant/Liner/Propellant  Bondedjpecimen/-  -  AdA-  4U  f 


In  this  sub-task,  numerical  modeling  analysis,  based  on  an  elastic-plastic  analysis,  of  the 
propellant/liner/propellant  bonded  specimen  will  be  conducted.  The  information  obtained  from  the 
numerical  analysis  will  be  used  to  explain  the  phenomena  observed  from  experiments.  At  the  time  of 
this  writing,  a  finite  element  model  was  developed  and  a  detailed  analysis  was  started. 
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Fig.  2  Comparison  of  Stress-Strain  Curves 
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Fig^  Photoelastic  Fringe  Patterns 
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